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THE DYNAMIC COUPLING BETWEEN 
F2-REGION AND PROTONOSPHERE 

H. G. Mayr* 

ABSTRACT 

The time dependent continuity equations for H + and 0 + are solved 
by means of a perturbation method. A solution for the diurnal varia- 
tion of these ion species is presented to illustrate the dynamic effects 
of the coupling between F2-region and protonosphere. The variability 
of the ion composition transition level is discussed in view of the ob- 
served ion trough at high latitudes , and from this evidence is found for 
a continuous depletion of protons there. 


*NAS-NRC Resident Research Associate, now with NASA. 
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THE DYNAMIC COUPLING BETWEEN 
F2— REGION AND PROTONOSPHERE 


INTRODUCTION 

Hanson and Ortenburger (1961) were first to discuss the coupling, or 
more appropriate the lack of coupling between F2- region (where oxygen ions, 
0 + , are abundant) and protonosphere (where hydrogen ions, H + , are dominant). 
They recognized the very significance of the diffusion barrier, H + encounters in 
collisions with 0 + when diffusing into and out of the protonosphere. This, and 
considering the fact that only a limited amount of protons is supplied through 
the charge exchange reactions 


H + ^ 0 ^ H + 0 + (1) 

(which are the most significant chemical processes for H + ), Hanson and Orten- 
burger predicted that the characteristic time for refilling and depleting the 
protonosphere is in the order of a day and therefore one should expect insignifi- 
cant diurnial variability of the ion content in the protonosphere. 

As consequence of this inhibited and slow process of ionization exchange 
between F2-region and protonosphere, the assumption of steady state in describ- 
ing these two regions (which would demand a short response time) is not valid. 
Their time dependence has to be considered as it affects the spatial structure 
of the ions. 

Subject of this paper is to present a solution of the time dependent con- 
tinuity equations for H + and 0 + in the topside ionosphere. A comparison between 
the steady state and time dependent solutions will allow to set criteria for the 
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study of certain magnetospherie and ionospheric phenomena that are yet 
unresolved. ’ 

THEORY 

Hanson and Patterson (1964) and Geisler (1967) have demonstrated that 
0 + or H + distributions are only significantly affected by ion transport fluxes if 
these are in the vicinity of the critical fluxes. Critical fluxes are upper limits 
on ion fluxes that can be supported by chemical production and diffusion processes. 
For 0 + the critical flux is during day in the order of 5 x 10 8 /cm 2 see and thus 
significantly higher than that for H + which is close to the 10 8 /cm 2 sec (during 
minimum solar activity, Brinton et al. 1968). During nighttime however, the 
critical flux for H + may be comparable to that for 0 + . The consequence is that 
dynamic variations in the protonosphere induce proton fluxes that, through a 
chain of chemical reactions also induce Q + fluxes and they may affect the density 
of the oxygen ions. For this reason it may not be justified to treat the time 
dependent variation of the F2- region isolated without considering the proton- 
osphere especially during nighttime. 

Extensive theoretical work was done on the dynamic of the F2- region. 

Ruster (1967) solved the time dependent equations of motion and particle con- 
tinuity to describe temporal fluctuations in the F2-region. Kritt (1967) solved 
the time dependent continuity equations and the quasi- stationary motion equations 
analytically to describe the diurnial variations of the 0 + concentration. Theo- 
retical models of the F2- region have the advantage that they are linear because 
the ion concentration is several orders of magnitudes lower than the density of 
the neutral atmosphere which thus determines the diffusion coefficient. In the 
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topside ionosphere, where the transition occurs into the protonosphere, the neu- 
tral atmosphere plays a minor role and the diffusion coefficients for 0 + and H + 
depend primarily on the concentrations of these two ion species. For this 
reason we deal with a nonlinear diffusion problem which is mathematically 
difficult to solve. 

In the following it is assumed that the diurnal variation in the F2- region 
is known at 300 km and it is asked for the diurnal variation of 0 + and H + above 
this altitude. The continuity equations for the two ion constituents are 


<?[ 0 + 
r) t 


P 0 - -f R[HJ f.0 + ] + RlH + ] LOl 


B 


<?s 


0 lH + ] 
dt 


9 _ 

8 


R[H] [0 + ] - R[H + ] [0] - 



(0) 


if diffusion along field lines (with distance s) is assumed to be dominant. It is 
R charge exchange coefficient (10 -1 1 T 1/2 cm 3 /sec) 
t time 

V macroscopic transport velocity parallel to the magnetic field 
The function B is proportional to the magnetic field 


B ^ ~ (1 ■+ 3 sin 2 A.) 1/2 

r J 

with 

r geocentric distance 
A geomagnetic latitude 
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and accounts for the divergence of the magnetic field tubes. The first terms in 
(2) and (3) express the photoionization rates (P 0 a [0], P H a [HJ ) the 
next two terms the charge exchange process (1) and the last terms stand for the 
transport process. 

Adding Equations (2) and (3) leads after integration to 


[H + ] V . + [0 + ] V , 

H 0 + 



d [o + l (Hh + A 
at " ot ) " c 


(4) 


s e is the field line distance to the equator. In Equation 4 some of the following 
boundary conditions 




d [Q + ] 
ds 


d [H + ] 
d s 


0 , 


for 


s = s e (5) 


were applied which lead to solutions that are symmetrical with respect to the 
equatorial plane and thus apply to equinox. 

The linearized equations of motion are for ions 


<9V 


m + dt + 6 + 

O ° L H 


+ [H + ] [0 + ] N - V + ) + 9 [0 + ] N V 

0 + ^ 0 + H + / 0 + n n o 


kT ~ ^° s ] + [0 + ] m + f + e[0 + ] E ^ 
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m . LH + ] 

H + 


dV + 

H 

d * H + 0 + 


[H + ] L0 + ] 




+ o , Lh + ] n V . 

H + n n H + 


= - kT ^ * [H + ] m + f + e t H + 1 E (7) 

and for electrons 

( 8 ) 

- kT -jj - cN e E ' 0 W 

where 

N n neutral density 

T electron temperature, which is assumed to be equal to the ion tempera- 
ture and to be constant with altitude, 
f acceleration parallel to the magnetic field due to gravitational and 
centripetal forces 

E electric polarization field parallel to the magnetic field 
k Boltzmann constant 
m mass 

e electric charge 
electron density 

The coefficients 6 are related to the diffusion coefficients D accord- 

x y * y 

ing to 


" ( CX] + [Y]) D xy ’ 

Estimations reveal that the time derivatives in Equations 0 and 7 are orders of 
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magnitudes smaller than other terms, if diurnal variations are considered. 
Therefore the motion equations can be treated as quasi- stationary. 

We eliminate the electric field by substituting Equation 8 into ((>) which 
yields 


e 


+ [H + ] [0 + ] (V - V + ) + 0 + [ 0 +] N v 4 

O \ O H / O n n 0 + 


= - kT 


<9[0 + ] [0 + ] dN 


ds 


kT 


+ [0+J m n+ f (9) 

e O 


Introducing (4) into (9) results in 


[0 + ] V 

0 


(0 . , N + 0 . N \ 

[ H + 0 + e 0 + n n ) 


'kT 


d [o + ] 

d s 


+ kT 


[o^j aw. 

N„ d s 


+ [O 


("V 


f + 0 


h + o h 


r do) 


In Equation (2) we substitute V + [0 + ] through (10) and integrate it by applying 
condition (5), This leads to the integral differential equation 


kT 


d [o + ] Lo + J * N . 


ds 


+ kT 


N e ds ~ m Q 


, f [o + ] = - e cLo + ] + (0 n 

+ H + 0 + \ H + 0 4 ’■ 


+ Vn N ") B J ^( P o " 8 KML° + ] + RlH + j [0 ] ) 


(11) 
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and correspondingly for H + to 


d [H + ] 


+ kT 


[H + l * N « 


N 


ds 


m + f[H + ] 


0 + + C[H + J + (Q + + N 

HO \ H O e 



(?[H + 1 

dt 


r[h] lO + ] - R[H + ] to] 


( 12 ) 


Considering that N e - [0 + ] + [H + ], the Equations (11) and (12) describe the 
temporal and spatial distributions (along field lines) for H + , 0 + and N ( , . 

As boundary condition for integrating (11) and (12) we assume that at 300 
km the temporal variation of 0 + is known and has the form 


[ 0 + ] 


= [o + l 0 + to + ] 



(13) 


where t is the period of one day. H + is assumed to be in chemical equilibrium 
through charge exchange and therefore is related to 0 + , 0 and H by 


Lh + ] 


9 LO + ] [H] 

8 to! 


up to 400 km. 0 and H are assumed to be in diffusive equilibrium whereby their 
values vary at any altitude with the same period and phase as does 0 + (Equation 
13). Also for T and T n the same form of temporal variation is adopted. 


SOLUTION PROCEDURE 

The Equations (11) and (12) are solved by means of an iterative approxima- 
tion in which the terms on the right hand side are treated as perturbations. For 
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stationary conditions (d/d t = 0)this method was employed in (Mayr et hi. 

(1967)), where it provided solutions accurate to 1% with computation times less 
than 10 sec on the IBM 360-75. We have therefore extended this method to solve 
the time dependent problem. This was accomplished by treating the time deriva- 
tives in (11) and (12) as perturbations which were applied iteratively in the fol- 
lowing scheme. 

With the k-th solution the time derivatives are assumed to be 



where e k is a parameter that varies between 


For the initial solution it is assumed that e = 0 which is equivalent of assum- 
ing steady state. In subsequent solutions e k is slowly increased which means 
that we allow the time derivatives to become slowly effective. This assures the 
applicability of the perturbation method which demands small effects from the 
"perturbing" time gradient. A successive application of this procedure leads 
then to a convergence on the solution of the time dependent continuity equations 
when € k approaches 1. 
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RESULTS AND DISCUSSION 

In the following, solutions are presented and discussed that were computed 
along a midlatitude field line that crosses 1000 km at 35°. Table 1 shows the 
time average (with subscript 0 ) and the amplitude (with subscript 1 ) for the input 
parameters. In Figures 1 through 3 the time dependent solutions for 0 + , H + and 
N g are shown in solid lines and for comparison the quasi- stationary solutions 
(d/d t = 0 ) are shown in dashed lines. 


Table 1 

List of Input Parameters. 


lO + ] 300 (cm" 3 ) 

[0 + ] 0 = 2.0 x. 10 s 

[0 + ] j = 5.0 x 10 4 

tO] soo (cm' 3 ) 

[O ] 0 * 1.5 x 10 7 

[OJ t = 1.0 x 10 7 

[H] 500 (cm -3 ) 

[H] 0 - 2.0 x 10 s 

tn] i - 2.0 x io 4 

T (°K) 

T 0 = 2.0 x 10 3 

Tj * 5.0 x 10 2 

T n (°K) 

T n0 1-0 x 10 3 

T . 2.0 x 10 2 

n 1 


As evident, the differences between time dependent and stationary solutions 
are significant and thus the assumption of steady state is not appropriate in de- 
scribing the topside ionosphere. An interpretation of the various dynamic effects 
requires a thorough study of the complex interaction of the many physical quan- 
tities involved. As this is beyond the scope of this paper we shall confine our- 
selves to briefly discuss some of the major features. 

The 0 + concentration exhibits at all altitudes a maximum during noon. 

The amplitude of its diurnial variation increases strongly with increasing altitude 
which is primarily the effect of the temporal variation in the plasma tempera- 
ture. At lower altitudes, the differences between the stationary and time 
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dependent solutions are negligible. At 700 km and above differences become 
apparent, first only during night at higher altitudes even during day. Differ- 
ences of more than 30% and a negative phase shift are clearly evident at 1000 km. 

The H + density shows a significantly more complex behavior. At 300 km 
its peak is in phase with that of 0 + and its minimum stretches over a rather 
long period of 8 hours during night. This results from the ion chemistry that 
involves variations of 0 + , 0 and H. At 500 km, the temporal variation of H + 
shows an entirely different picture. A broad minimum is formed during day 
and a sharp maximum developes at night. The reason for this change can again 
be found in the ion chemistry which at this altitude involves significantly differ- 
ent [O] / [H] ratios. In addition it is evident from the difference between the 
stationary and time dependent solutions that downward fluxes of II + enhance the 
proton concentration at night by almost a factor of two. At 700 km the temporal 
variations of H + result from the interaction of several processes. There the 
transition occurs from the F2-region into the protonosphere which means 
for H + that transport processes start to overrule chemical processes. 

This is apparent in the large differences between the stationary and time de- 
pendent solutions which show that during day upward fluxes decrease the proton 
density while at night downward fluxes enhance the H + concentration. These 
features prevail up to 2500 km where the H + density is decreased during day 
and enhanced at night, such, that the amplitude of the temporal variation de- 
creases by a factor of two and the phase shifts by almost three hours. 

The dynamic effects are clearest apparent in the electron density distribu- 
tion (Figure 3). It shows that with increasing altitude the amplitude decreases 
while density minima and maxima shift in local time. These effects occur 
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during night between 500 and 1000 km and during day between 700 and 2500 km, 
the reason being that the 0 + -H + transition region is at night at a lower altitude 
than it is during day. 

In Figure 4 the 0 + -H + transition level is shown for the time dependent 
solution in solid line and for the steady state solution in dashed line. The dif- 
ference is very pronounced, showing that dynamic effects enhance the transition 
level during day by 150 km and decrease it at night by 100 1cm. Thus the ampli- 
tude of the variation of the transition level increases as the result of dynamic 
effects. This feature may provide a simple qualitative criterion to test certain 
ionospheric and magnetospheric phenomena. 

Ion composition measurements (Taylor et al. 1968) have shown that at high 
latitudes (> 50°) the mean ion mass at 1000 km is significantly higher than the 
value one should expect on theoretical grounds when assuming the ion composi- 
tion to be in steady state (Mayr et al. 1967). This discrepancy can be interpreted 
as due to upward fluxes of protons that increase the mean ion mass and the 
0 + -H + transition level. As shown previously, during some period of the day such 
fluxes could result from the dynamic coupling between iononosphere and proton- 
osphere. But, if this is the case then downward fluxes should prevail during 
some other period of the day and these fluxes would decrease transition level 
and mean ion mass below their steady state values. Taylor's mean ion mass 
observations were made during dawn and dusk periods and they do not show 
significant differences thus suggesting that the previously discussed dynamic 
effects may not be dominant. Unless a more detailed experimental coverage 
of the diurnal variation in the ion composition reveals that between dusk and 
dawn significant variations of m + and transition level exist at high latitudes, 
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one has to postulate a continuous upward flow of protons. This would indicate 
that plasma is leaking across the magnetic field or escapes in open field tubes 
into the interplanetary space, mechanisms that were employed to explain the 
plasmapause and its relation to the ion composition at higher latitudes (Mayr 
1968 ). 
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FIGURE CAPTIONS 


Figure 1 Computed temporal and spatial distribution of 0 + . In solid lines the 
time dependent solution, in dashed lines the quasi- stationary solution 
( d/dt = 0). 

Figure 2 Computed temporal and spatial distribution of H + . In solid lines the 
time dependent solution, in dashed lines the quasi- stationary solution. 

Figure 3 Computed temporal and spatial distribution of N e . In solid lines the 
time dependent solution, in dashed lines the quasi- stationary solution. 

Figure 4 Computed temporal variation of the 0 + -H + transition level. In solid 
line the time dependent solution, in dashed line the qua si-stationary 
solution. 
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Figure I. Computed temporal and spatial distribution of 0 + , In solid lines the time dependent 
solution, in dashed lines the quasi-stationqry solution (rl rh 0 ). 
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Figure 4. Computed temporal variation of the 0 + -H + transition level. In solid line the time 
dependent solution, in dashed line the quasi-stationary solution. 
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